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1. Introduction 

Experimental setups in ultracold atomic gases allowing easy optical access to the 
sample and the well established imaging techniques have resulted in a wide range of 
precise and well understood spectroscopies. Radio-frequency [TJ [^J [31 SI El El [3 IS] (RF), 
Raman [9] [10], Bragg [TTJ [12l E2 03] , and lattice modulation spectroscopies [I51H61IT7] 
have plenty in common and they can all be described effectively as an absorption of a 
single photon of momentum k and frequency uj. 

The agreement between spectroscopic experiments and theories in ultracold gases 
has often been at most qualitative despite the various theoretical approaches such as 
mean-field [UmilBISll^llHIMlHSlllS, t- matrix [37], or other diagrammatic 
approaches [28] [29] [30] [31] . Strange peaks and features (such as anomalously broad 
spectral peaks) in the experiments have been attributed to the lack of coherence [14] . 
defects in a lattice [15] , or to center-of-mass motion of pairs [8] , suggesting that those 
features are caused by imperfections in the experimental setups. These anomalies 
have been all but neglected in theoretical descriptions and, instead, a lot of effort has 
been put into improving the underlying many-body theories by including thermal and 
quantum fluctuations [32, 33, 34] in order to obtain quantitative agreement with the 
well understood features in the spectra. Here we show that those neglected anomalies 
seen in a wide range of different experiments have a common source. 

While the breakdown of a mean-field theory in strongly interacting systems is 
fully anticipated, the treatment of beyond mean-field effects is exceedingly difficult. 
As a result most theoretical studies have relied on the mean-field approximation, 
partly because of the apparently good agreement between theory and experiments. 
However, as we will soon point out, various experiments performed with strongly 
interacting systems have features that cannot be explained by mean-field theories. 
Here we propose an efficient diagrammatic expansion for the mean-field theories to 
describe an important class of coherent processes. 

2. Linear response spectroscopy 

We begin by deriving the linear response theory for a weak perturbation V. While 
the calculation is a simple textbook derivation, it will help to illuminate the problem 
encountered with the mean-field linear response theory. 

Let's assume that the system is initially (at time t = 0) in a ground state of a 
(many-body) Hamiltonian Hq, denoted by \i/)(0)). Switching on perturbation V(t) 
yields the full time-dependent Hamiltonian H(t) = H n + V(t), and the system evolves 
according to the time-dependent Schrodinger equation 

ihj t m))=H{t)\m)- 

This can be solved to yield the state at time t as 

\m) = e- iH ^ H s(t)\m), 

where the scattering S-matrix is 

s(t) =1+4 w) + ( 4i f dt> f ^ vitfw) + . . . , 

in Jo \ lh J Jo Jo 

where V\{t) = e lH ot/fry^ e -iHot/fr j s ^ ne perturbation expressed in the interaction 
picture. This formulation of S-matrix readily offers easy access to the perturbative 
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description. Cutting the series expansion of the S-matrix yields perturbation theories 
of different orders. For our purposes, it is enough to consider the first three terms 
of the above expansion, yielding the approximate wavefunction in the second-order 
perturbation theory 

2 ,4 r t> 



\m) 



-iHot/h 



1 



i / dt'V^+f^-] [ dt> [ dt" 
ihj \inj J Q J 



1^(0)) . (i) 



Ultimately, we are interested in the expectation value of some observable O. 
Keeping terms only up to the second order in the perturbation V, we obtain 



(o(t)) = m)\o\m) = (o) 



i 



dt 1 / a? (ViWawvitf')) 



W 1 dt ' 



df KOi(i)vKfw(t")> + www))] , 



(2) 



where the expectation values on the right hand side are calculated for the initial state 
\ip(0))- Should we be interested in the rate of change of (0(t)) (as is the case when 
measuring current), one can differentiate Eq. ^ with respect to time t and obtain 

j t {0{t)) =^jyt'^[{V I (t)0 I (t)V I (t')) (O^WWiit'))} , (3) 

where 5ft is the real-part and we have assumed here that V(t) and O are Hermitian 
operators. In order to continue this discussion on linear response, we need to specify 
the perturbation V and the observable O. 



2.1. Radio-frequency spectroscopy 

We will first consider the case of radio-frequency spectroscopy. A radio-frequency 
field of frequency ui and amplitude £1 drives atoms from hyperfine state |2) to initially 
unoccupied hyperfine state |3). As the wavelength of the radio- frequency field is much 
longer than any relevant length scale in the system, the radio-frequency photons do 
not affect the momentum of the atoms but only couple to the hyperfine state degrees 
of freedom. The perturbation in the rotating wave approximation is then given by 
V(t) — £le l5t / h ^2 k Cg fe C2fe + H.c, where 5 is the detuning of the rf-field frequency away 
from the Zeeman splitting of hyperfine states \2) and |3), S = fku + E^ f — E? { , where 
E?j is the hyperfine energy of state \a). The observable of interest is the number of 
atoms in hyperfine state |3): O — N 3 = J2k c 3k c 3k- Since the hyperfine state |3) is 
initially empty, -ZV3 |i/j(0)} = 0, the formula for the current simplifies into the form: 

— (N 3 (t)) = 4 I dt' me lHot/h V(t)N 3 e~ lHo ^ t,)/h V(t')e- lHot ' /h )JA) 
dt h 2 Jo 

where we have written the operators in the Schrodinger picture. Operating by .ZV3 
within the expectation value on the right hand side of Eq. (j4} will necessarily give 1 
since the single perturbation V{t') has transferred only one atom into the hyperfine 
state 1 3) and we assume that A3 commutes with Ho- Thus we have 

±(N 3 (t)) = ^^^^5ft e -^(-OA (c t pC3pe -^( t -OA4 C2fc)j 
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where we have assumed that the initial state |^(0)) is an eigenstate of H with 
eigenenergy equal to zero (a nonzero energy would provide only phase factor). This 
can be further simplified to yield 

j t (N 3 (t)) = ^5> f\t'e-^^(cl pC3p e-^'/% k c 2k ). 

kp 

The expectation value on the right-hand side is in principle a two-particle (or particle- 
hole) Green's function propagating within a many-body sea of fermions. However, 
in the special case where the hyperfine state |3) is noninteracting the two-particle 
propagator separates into a product of one-particle and one-hole propagators 

^W)) = ^^^dt'e- iSt '^G 2 (k,-t')G 3 (k,t'), (5) 

where —iG2(k,—t') — {c\ k e~ lHat l n c 2 k) is the propagator of a hole in single-particle 
state |2, fc) and iG 3 {k,t') = (c 3 ke^ lHot ^ hc \k) IS ^ ne propagator of an atom in state 
|3, k). Being a noninteracting Green's function (and since hyperfine state |3) is initially 
empty), the particle propagator G 3 (k,t') is simply —ie~ ltkt / h , where is the energy 
of an atom with momentum k and we end up with 

j t (N 3 (t)) = ^^jy t 'e- i ^ t '/' l iG 2 (k,-t') 

In the limit of a very long pulse t — > oo (but keeping the perturbation weak tft/h <C 1), 
the integral yields the Fourier transform of the real-time hole propagator. 

|TO)) = ^]>>G 2 (fc,e fe + *), 

k 

where 3 is the imaginary part. Thus the radio- frequency spectrum probes the spectral 
function of the hole A 2 (k,e) = 23G 2 (fc,e) which is ultimately determined by the 
operator e ~ lHat / h . 

While the simple BCS mean-field treatment has proven to be able to describe 
spectroscopies of even strongly interacting gases rather well in many cases, quantitative 
agreement between theory and the experiments has turned out to be very hard to 
achieve. Therefore, it is worth looking more closely what information has been lost 
in the BCS mean-field approximation. In particular, we are looking for non-BCS 
processes which are important for the hole propagator in order to obtain quantitative 
agreement with the observed spectra. 

2.2. BCS mean-field approximation 

We want to describe radio-frequency spectroscopy of a two-component Fermi gas. 
Since the rf-field couples one of the two hyperfine states into a third excited state, 
the full system is described by a three-component Hamiltonian but in which the third 
state 1 3) is initially empty. The Hamiltonian of the three-component system is 

Ho = ^2 e k c ka Ck ° + 2 9 arT ' C k,a C -k+q,a' C -P+Q^' C P,< :r ^ 
ka kpqa^a' 

where €k is the single-particle kinetic energy, a £ {1,2,3} and g aa i is the interaction 
strength between the atoms in hyperfine states \a) and |er') in contact interaction 
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potential approximation. The standard approach uses the two-body scattering T- 
matrix with g aa > = 47r ^' T ' T ' y, where a CTCT / is the scattering length and V is the 
volume. Initially we will assume that only 1712 is nonzero but later on we will consider 
also the general case in which all interactions can be nonzero. 

The BCS mean-held approximation has proven to be very simple but still very 
effective in describing two-component Fermi gases even in the very strongly interacting 
regime. The approximation amounts to first considering only opposite momentum 
scatterings q — in the Hamiltonian Hq and then replacing the quartic interaction in 
Hq by quadratic mean-field interaction, resulting in BCS mean-field Hamiltonian 

#o BCS = E £ * c L<*«t + A u J2 4,i c -*,2 + H - c -> 

ka k 

where the pairing field A12 = 1712 X) P ( c -p,2Cp.i)- Here we have neglected the Hartree 
energy shift but it needs to be included when studying for example gases in optical 
lattices in order to obtain correct bound state energies. 

The BCS mean-field theory can be solved exactly by diagonalizing the BCS 
Hamiltonian using the Bogoliubov transformation. This yields the Hamiltonian in 
the quasiparticle basis 

ka 

where C is a constant and the operators jk.a (7I a ) annihilate (create) quasiparticle 



excitation of momentum k and energy Ek = y (e/- — fj,) + Af 2 with /i being the 
chemical potential of the atoms in hyperfine states |1) and |2), and e k is the single- 
particle excitation energy. The quasiparticle 'spin' a € {t>4-} does not refer to the 
hyperfine level of an atom but simply to the fact that there are two branches of 
quasiparticle excitations. The hole spectral function can be shown to be 

e-E k + ir) e + E k +ir) 

where the Bogoliubov coefficients u\ = ^(1+ and = 1 — v%. We have assumed 

that the system is in equilibrium and the quasiparticle distributions are given by the 
Fermi-Dirac distribution np (e) = 1/(1 + e^), (3 = l/(k B T), k B is the Boltzmann 
constant, and T is the temperature. We consider only zero temperature T = results 
in this manuscript. 

At zero temperature, the BCS hole spectral function (|6]) consists of a single 
sharp peak (second peak appears at finite temperatures), describing a quasiparticle 
excitation with infinite lifetime. This property of infinite lifetime can be easily seen 
from the form of the BCS Hamiltonian H® cs . Creating a hole in state \k,2) will 
necessarily render any atom in state | — k, 1) effectively noninteracting, since the 
Hamiltonian H BCS describes only scatterings between atoms with opposite momenta. 

In order to consider quasiparticle lifetime effects, one would need to go beyond 
the BCS approximation and include residual interactions between the quasiparticles 

H r Q cs =g 12 c ii c -k+ q ,2 c -p+i^c p ,i (7) 

kp.q^O 

This approach has been studied in the context of nuclear physics [35] and high- 
temperature superconductors [3 6) . It has also been discussed in the present context in 
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Ref. [37] using a fully self-consistent calculation of the many-body T-matrix. The self- 
consistent nature allowed the inclusion of Anderson-Bogoliubov (AB) phonon effects, 
and these were shown to result in broadening of spectral features. The theory was, 
however, highly involved, and the present theory can be seen as an attempt to include 
important quasiparticle lifetime effects in a simpler way, namely using first-order 
perturbation theory. 

While Eq. (J7J for residual interactions is enough for our purposes, it is very 
instructive to briefly consider its presentation in the BCS quasiparticle basis. Writing 
the residual interaction Hq cs — -^40+^31 + ^22 in the quasiparticle basis yields [351157] 

H40 = .912 ^2 Vk+ q V k U v U p+ <f{k,-Xl-k-q,l1-p,i"(v+qA + ( 8 ) 
kp,q^0 

H 3 1 = 312 ^ (v p V p+g V k U k - q - U p U p+q UkVk- q )ll a lk^q,al- P ,ilp+q,t + H.C., (9) 
kp,q^0,a 

and 

#22 = 512 [(uq-kUkU p Uq-. p + v q - k v k v p v q - p ) 7g-fc, t 7fc,4.7p,4-7g-p,t 

+ (u k U p - q V p V k -q + Vk.Vp-qUpUk-q) ll^l^ p ^lq- P ,ilk-q,\ 



+Uk+qU p+ qV k V p jl +qt -fl t J p+ q^k,t 
+ Vk + qV p+ qU k U p -fl + qi -flp p+ q illK l 



(10) 



Clearly the number of quasiparticles does not need to be conserved as a single 
excitation can be split into several excitations through scattering and vice versa. 
Moreover, since quasiparticles have minimum energy cost of A, transitions to states 
with different numbers of quasiparticles are gapped. 

However, these quasiparticle interactions provide also coupling to the collective 
Anderson-Bogoliubov (AB) mode which is a phonon-like excitation and thus gapless. 
The description of this mode requires nonpcrturbative treatment, for example random 
phase approximation (38] [39] , and while the AB phonon arises already from the 
mean-field theory, the actual coupling between the quasiparticle excitations and these 
collective modes is provided by the residual interactions Hq 6S [57] . What our theory 
does is to include approximatively the coupling between the quasiparticle excitation 
and the phonon modes. We will later consider also settings in which the couplings to 
different numbers of quasiparticle excitations will be important and such processes 
may actually arise naturally from the perturbative treatment. Furthermore, we 
will also consider systems at temperatures above the critical temperature where the 
supcrfluid order parameter vanishes. In such case, the BCS model yields an ideal finite 
temperature Fermi gas without pairing correlations, binding energies, or Anderson- 
Bogoliubov phonons. In that limit, particle and the quasiparticle become the same. 



3. Beyond BCS mean-field approximation 

Our goal is an effective theory which can include the effect of interactions between 
quasiparticles approximatively and which can be easily incorporated into the standard 
mean-field BCS theory. As our starting point we take the residual interactions 
neglected in the BCS mean-field theory and treat these as a perturbation in 
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the hole propagator. Writing the time-evolution operator e lH ° f in Eq. ([T} in terms 
of the scattering S-matrix as 

e -iH t/H = e -^ CS */^ rcs(i)j 

and expanding the S-matrix S ICS {t) = 1 + i /„' dt' H*f(t') + where JZJf (t') = 
e iH™ cs t jjres e -iH™ cs t ^ we o^ain the first-order correction (in terms of i?o° s ) 

= (Jj^ ftff dt" e -^ cs (t-t')/n 

~H r es e- iH °° S (*'-*")/ + ^') e -i"o BCS (f'-n/^es 

e-^? 08 *"/' 1 ^^)) (11) 

Here we have included only the first order term (with respect to rf-coupling V) from 
Eq. ((T|) since the zero and second order terms do not contribute to the current in 
Eq. (U). The above equation is a result of applying perturbation theory to both 
residual interaction Hq bs and the rf- field V(t). 

At this stage, we notice the similarity with two-photon processes in quantum 
optics, such as the two-photon Bragg scattering or two-photon Raman transitions 
which also incorporate correlated couplings with two different fields. In the present 
case the two fields are the rf- field and the atoms in the hyperfine state |1). The 
main difference, however, is that in our model the atom-atom interaction has many 
different scattering channels, including very long wavelength (low momentum exchange 
q) transitions. Nonetheless, the similarity is encouraging - we are here dealing with 
physics that has manifested in other systems as well. 

The first term in Eq. (fTTj) can be depicted pictorially as a diagram in Fig. [1] (the 
second term corresponds to diagram in Fig. [3] as will be discussed later). The diagram 
has been drawn for a more general spectroscopy than rf-spectroscopy by including 
also the photon momentum q^. Below we will assume that q^ = for rf-spectroscopy, 
but later on in the discussion regarding Bragg spectroscopy the momentum of the 
photon can be significant. The diagram describes a path for the wavefunction and one 
will need to calculate the squared norm of the amplitude in order to find expectation 
values, such as a Green's function. As in the discussion regarding linear response 
above, a first-order perturbation for the wavefunction corresponds to a second-order 
perturbation for the propagator. 

In the case of rf-spectroscopy, the diagram describes a correlated process in which 
a rf-photon creates an atom in state |3, k) and a hole in |2, k). The hole then scatters 
with a hole in state \l,p) into a new set of holes |2, k + q), \l,p — q). Another way 
to describe the diagram is that an atom in state |2, k) absorbs the photon and is 
transferred into the state |3, fc). After the rf-photon absorption, another pair of atoms, 
originally in states |2, k+q) and \l,p—q), scatter into states |2, k), \l,p), filling the hole 
created by the RF-transition. This description can be better illustrated with diagram 
in Fig. [2] While the latter description sounds like the processes (rf-photon absorption 
and atom-atom scattering) are disconnected, the first interpretation shows that they 
are indeed connected through hole degrees of freedom. However, below we will use the 
latter description as it generalizes more easily. The corresponding diagrams for holes 
can be easily formulated as the two representations are completely interchangeable. 

In the following, we will assume that the propagator lines for atoms in hyperfine 
state |1) are described by bare Green's functions and only the atoms (or holes) in 
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912 



< 



1,P 



Figure 1. The lowest order photon absorption Hlu, and quasiparticle scattering 
912 process. The absorption of the photon by an atom in the ground state of |2) 
produces an atom in hyperfine state |3) and one hole (backward propagator) 
excitation in |2). The hole excitation in |2,fc) then scatters with a hole in 
resulting in holes at states \2,k + q) and \X,p — q). As is clear from this 
description, the photon absorption and the atom-atom scattering processes are 
connected through the hole degrees of freedom. If the photon has finite momentum 
9l such as in Bragg spectroscopy, the atom will also feel a momentum change 
worth ql- We approximate that |l)-atoms are described by bare propagators and 
hence the atom-atom scattering process does not produce additional excitations 
experiencing mean- field energy shifts (such as BCS pairing gap A). The process 
describes qualitatively the coupling of the hole excitation to a gapless phonon 
mode. 




3, k + g L 



2,k + q 



2,k 





512 1 






^ — 


> 



Figure 2. An alternative way of seeing the process in Fig.[T] The absorption of 
the photon by an atom in the ground state of |2) produces a particle in hyperfine 
state |3). After the photon absorption, another pair of atoms in states |2, k + q) 
and |l,p — q) scatter into states |2, k), filling the hole created by the photon 

absorption. 



hyperfine state \2) are described by dressed propagators. If we look at the residual 
interactions H}f s in the quasiparticle basis, it is clear that in the first order of scattering 
the single quasiparticle created by the RF-field couples only to a state with three 
quasiparticles (through Eq. ((9|)) and not to the gapless Anderson-Bogoliubov phonon 
with a single quasiparticle. This problem is generic to a perturbative approach, 
as the description of the collective Anderson-Bogoliubov phonon would require self- 
consistent, or nonpcrturbative, treatment. The assumption invoked here allows us to 
capture elements from the Anderson-Bogoliubov phonon since the scattering does not 
introduce extra quasiparticles. Indeed, in the zero-momentum scattering asymptote 
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q — > 0, the energy shift due to the scattering becomes Hcq, where c = — . Close to the 
Fermi surface, p~k F , this yields the speed of sound as expected for the AB phonon. 

Thus by using bare propagators for atoms in state |1), the diagram [2] yields an 
approximative description for the coupling to the Anderson-Bogoliubov mode. To 
elaborate, by allowing the single quasiparticle excitation to scatter into lower energy 
states, the hole wavefunction will spread over several momentum states leading into the 
decay of the excitation. However, since the process is correlated with the rf-photon 
absorption it is possible for the quasiparticle excitation to scatter also to a higher 
energy state if the rf-photon has excess energy. Simply put, this means that the excess 
energy of the rf-photon is transferred into extra kinetic energy of atoms. However, 
since the rf-transition conserves momentum, this change in kinetic energy is provided 
by the atom-atom scattering. Similarly, given extra energy from the probing photon, 
a quasiparticle at the bottom of the quasiparticle energy band (corresponding to a 
hole at the Fermi surface) can decay by scattering to higher energy states, acquiring 
a finite lifetime. This might be connected to the observation of a finite linewidth of 
such quasiparticles in Ref. [37] . 

The diagram [2] can be added to the standard linear response calculation of 
the radio-frequency spectroscopy, thus generalizing the linear response calculation to 
include couplings with gapless phonon modes. As discussed already in Ref. [37], such 
coupling broadens the sharp quasiparticle peak in the hole spectral function. We will 
discuss the implications later on in the section [5] 

Including the diagram in Fig. [2] yields the rf-transition (for zero photon 
momentum <7l = 0) 

dt J\t'e^ +5E ^ {g 12 cl p 4 >k c 2 , k+q c 1>p ^ e -«^(*-*') (n C \ k c %k ) |^(0)) e -" t ', 

where SEi = e k — [i + El — S is the energy change due to rf-photon absorption 
(describing the energy of a single particle in state \k, 3) and hole in |2,fc)), SE 2 = 
€p — e p _ q + El +q — El is the energy change due to the scattering (corresponding to the 
kinetic energy difference of the atom in hyperfine state |1) and the change in the hole 
excitation energy due to the scattering of the hole from state |2, k) to state \2,k + q)), 
and the factor e~ nt has been added to guarantee convergence. The factor 77 plays 
the role of linewidth of the rf-field and it can be interpreted as a switching off of the 
rf-ficld. Here wc have included the hyperfine state indices in the single-particle and 
quasiparticle energies for clarity. The time integrals will yield the energy conservation 
but in order to calculate the actual transition probability amplitude, one will need to 
calculate also the norm of the state 

= C l,p C L C 2,fc+9 c l,P~94,fe C 2^IV'(0)), 

which depends on the details of the initial state |^(0)). The transition probability can 
be obtained by calculating the norm of this transition. However, as will be seen later, 
there are also other diagrams describing the same transition and these will need to be 
summed coherently before calculating the probabilities. 
The time integrals can be easily evaluated, yielding 

p ho i 1 / \ _ n mA, P 4,k c ^k+ q cx,p- q cl k c 2 ,k\^)) 

Mc.p.qlW 



5Ei + irj' 6 Ei + 5E 2 + in 

ft 9l2cl p cl k C2,k+ q Cl,p- q cl k C 2 ,k\^(0)) 

4~tJ< + Et-5 + ir)> el + el + E 2 k+q — fi — e p _, - 5 + ir, : 



(12) 
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Figure 3. The other first-order diagram for the photon absorption hu,qi, and 
atom-atom scattering 312 process. Here the atom in hyperfine state |2) scatters, 
producing one particle excitation and one hole excitation. The particle excitation 
is then transferred to the hyperfine state |3) by the photon absorption. 



7 tun,q L 



2,k + q %i 3,k + q + q L 3,k + q L 

7^ ^ 



913 



hp-q 



Figure 4. The first-order diagram for the rf-photon absorption huj, q-^ and for 
atom-atom scattering #13 process. Here again the rf-photon creates a particle in 
hyperfine state |3) and a hole excitation in |2). However, now the scattering takes 
place in the 13-channel. Assuming that initially the hyperfine state |3) is empty, 
this is the only first-order diagram contributing in this interaction channel. 



where we have added a lifetime broadening jf to the (virtual) intermediate state. We 
use the lifetime broadening rj = |gi2|iVi, where N\ is the number of particles in state 
|1), because the two-body scatterings are the dominant decay channel in a dilute gas. 

For narrow linewidth 77 <C rf , the transition probability amplitude has a single 
peak centered around 8 = (e| + e*) — (fi — E% + + ep_ g ). This corresponds to the 
energy change of a process in which one transfers one dressed ground state particle 
from state |2, k + q) and one bare particle from state \l,p — q) into two bare particles 
in states |3, k) and 1 1, p) . Notice that in the limit q — > the peak approaches the 
standard linear response result S = e| — fj, + E^, the asymptotic difference being linear 
in q as already discussed above. 



4. Sum rule considerations 



Based on very general arguments, the atom-atom interaction effects should vanish from 
the radio- frequency spectra when the interaction strengths between atoms in hyperfine 
states |1) and |2) (512) and between states |1) and |3) ((713) are equal (221 00]. Thus 
we expect that the diagram in Fig. [5] needs to be combined with other diagrams of 
the same order in photon coupling and atom-atom interactions. The other first-order 
diagrams are shown in Figs. [3] and [4] The contributions from these diagrams are 

„12 , A _ 912 ^ C 3,fc C 2,fc C l,p4,fc C 2,fc+gCl,p- g |V'(0)) 

+ £ p + ^k+q _ e p-q + l7 l e p + e k + ^k+q M e p-q ~ + tT) 
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13 , A n 9i3c\,pc{k c 3,k+ q Ci, P - q cl }k+q c 2 , k+q \Tp{0)) 

€ k+q M + ^fe+g - + "7 £ p + £ fe + ^/c+g M £ p-<? ° + iT l 

We refer to the diagrams [U |3J and 2] collectively by term resonant scattering and 
absorption (RSA) diagrams. 

The quasiparticle line entering the scattering vertex in diagram [3J describes a 
(mean-field) ground state pair (energy — E% + < 0) and the line leaving the scattering 
vertex is excited (E\ > 0). The role of the photon is now to transfer the particle 
excitation to the state |3). Since quasiparticle has already been excited in the 
atom-atom scattering process, the photon absorption itself does not create further 
excitations. This follows also from the observation that the total energy difference for 
all three diagrams is the same as they describe the same transition. The two diagrams 
shown in Figs. [2] and [4] do not create more than one quasiparticle excitation and hence 
the process in Fig. [3J cannot result in more than one excitation. 

The diagram [3J warrants even a closer examination. It describes a process in 
which the atoms scatter before the photon is absorbed. The first part of the diagram 
(the scattering) can be understood as a first-order correction to the mean-field ground 
state. Indeed, neglecting the photon line for a second, the state obtained from the 
first-order time-dependent perturbation theory, when the perturbation (in this case 
the atom-atom interaction) is switched on adiabatically, is the state one would have 
obtained from the first-order time- independent perturbation theory. The perturbative 
correction to the state can be understood as a fluctuation around the unperturbed 
state, and the photon in the diagram |3J couples to these fluctuations. While the 
diagram|3jcould be removed by including the first-order time- independent perturbative 
correction to the initial state \ip(0)) (that is, by including the fluctuations in the state 
|-0(O))), the diagram in Fig. [2] cannot. If we interpret the diagram [3] as a process 
in which fluctuations absorb the photon, the diagram [2] corresponds to a process in 
which the fluctuations are enhanced by the absorption of the photon. Clearly, the 
latter process cannot be included in the initial state |^(0)) as the photon field is 
absent at time t — 0. The present formulation as time-dependent perturbation theory 
allows an easy and consistent (in the sense of the sum rules discussed below) way to 
include the effect of fluctuations in the spectroscopic probes. 

4-1- Asymptotic high energy scaling 

The form of Eq. ([13")) has interesting implications for the asymptotic tail of the radio- 
frequency spectrum. In a uniform three-dimensional system the single- particle energy 
dispersion is €k = 2 • ^ n the high energy asymptote of the spectrum 8 3> Ep, 
pairing effects described by the gap A have little effect (allowing the approximation 

E\ = sj (ek — n) 2 + A 2 sa eu) and resonant processes are dominated by scatterings to 
high momentum states p, k 3> where kp is the Fermi momentum. Assuming that 
only low momentum states are initially occupied one has fc + ?~p-?~0 (hence 
k ~ — p and q is large). Now Eq. (fT3)) yields 

p 12-a S ymptote| _ gig " C 3,fc C 2,fc C L4,fc C 2,fc+gCl,p- g IV> (»)) 

In the large k limit, the probabilities for the states |2, k) and |l,p = —k) to be empty 
approach unity and likewise the occupation probabilities for states — q w 0) and 
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|2, k+q 0) approach unity. Calculating the transition probability I j^^"^ symptote | ^} | 2 
and integrating over k one obtains the transition probability for a single particle in 
state |2, k + q) to be transferred to hyperfine state |3, k) 

2 



pl2 — asymptote 



de- 



912 



2e + irf 2e — <5 + if) 



(16) 



where the ^fe factor comes from the density of states in a three-dimensional uniform 
system. For narrow linewidth rj -C 5, only energies e close to 5/2 contribute to the 
integral and it can be approximated as 



P 



12— asymptote 



512 



IT)' 



2e-5 



it) 



-, (gl2^) 7T 

(5 2 + rj' 2 ) rj 



asymptote 



(17) 



c-1.5 



Assuming 6 rj , we obtain the asymptotic high RF detuning tail P l 
It is important to notice that this result is independent of a possible 'quasiparticlc 
tail' in the occupation numbers W2,k- Diagrams [2] and 2] also produce high energy tails 
for the spectra. The asymptotic form of the latter diagram is 

2 



P 



13— asymptote 



de ■ 



n 



513 



(18) 



—5 + irf 2e — 5 + ir\ 

which yields the same asymptotic behaviour as Eq. (fl7j) but is proportional to 513 
instead of 512. Unlike the above two asymptotic forms, the high energy asymptote of 
diagram [5] requires occupation of high momentum states and it is given by 

2 



P 



hole— asymptote 



de ye 



(19) 



2e - 5 + irj 2e - 6 + in 
is the probability amplitude that there is an atom in hyperfine state 



For a BCS-like state this is given by the Bogoliubov coefficient 
A 2 ) yielding the same asymptotic high energy tail 6~ 



where v 2 
2 with energy e. 

v 2 e ~ 5 (l — e/ve 2 + A 2 ) yielding the same asymptotic high energy tail 5~ l ' b . Here 
we have neglected terms of the order v\ which is justified assuming that the occupation 
probabilities do not have too fat a tail. The most important point to notice, however, 
is that only Eq. ([T9")l depends on high energy excitations in the initial state IV'(O))- 

The 5~ 15 scaling law has been suggested already in various studies [TO 144) . 
where the scaling law comes from the high energy tail of the occupied atom states (the 
V2,k or U2, e tail). In present theory, no such many-body correlations are needed in the 
initial state - the power law arises directly from the RSA process. In particular, 
the RSA scaling law applies in the limit A —> and also to all momentum k modes 
separately, resulting in a (5~ 15 tail in the spectral function of every k mode. This will 
have implications in particular for the momentum resolved spectroscopy [8] that we 
will consider in Section UTTl 

All the RSA diagrams lead into the same final state \ip) and thus they interfere 
coherently. A closer examination of the orderings of the creation and annihilation 
operators shows that, because of fermionic anticommutation rules, Eq. (|14p has 
opposite sign to the Eqs. (| 12[) and (fl3|) . and hence the three contributions interfere 
destructively. This is easy to see by operating on Eq. (fT4)) with the identity operator 



1 



C 2,fc u 2./> 



+ C. 



:i,k c 2,k- 



In particular, for identical interactions (1712 — 513) and 
for equal dispersions {E\ — e|, for all k), the three RSA diagrams combine into a 
single sequential process. Such process does not allow energy exchange between the 
quasiparticle and the emitted phonon. Thus, in the symmetric case, our model does 
not provide any energy shift to the rf-spectrum. This is in agreement with sum rules. 
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The sum rule argument must hold for all orders of atom- atom scatterings, and 
hence it will hold also (or especially) in the nonperturbative treatment. Thus the 
picture from our model, in which the rf-photon absorption and the excitation of the 
phonon are correlated, should also hold for a real system in which the phonon is the 
Anderson-Bogoliubov phonon. 

If the final state |3) is interacting, also the BCS mean- field response changes 
and one needs to include the Aslamazov-Larkin (AL) and Maki-Thompson (MT) 
contributions. These resemble our beyond mean-field diagrams. However, these 
are strictly mean-field effects and do not contribute to the quasiparticle excitation 
lifetimes. The effect of these diagrams have been studied in the context of RF 
spectroscopy of atomic gases in Refs. [451131] . 



5. Comparison to experimental results: fermionic gases 

All spectroscopies used in the context of ultracold atom gases can be understood as an 
effective photon absorption and hence we expect diagrams similar to those in Fig. [2] to 
play a role. The main difference in various spectroscopic methods is the atom degree of 
freedom that is affected: hyperfme state (RF and Raman spectroscopies), momentum 
state (Bragg and Raman spectroscopies), or lattice band (lattice modulation and 
Bragg spectroscopies). Furthermore, the diagrams (with slight modifications) can 
be applied also to bosonic atoms. In the following we will analyze seven different 
experiments which we believe to show the effect of the RSA process. 



5.1. Radio-frequency spectroscopy of a Fermi gas 

An experiment performed in JILA in 2008 [HI HH] applied momentum resolved RF 
spectroscopy to probe the single-particle excitation spectrum of a two-component 
Fermi gas. The spectrum in Fig. 3c of Ref. [8] was shown to be broader than expected 
and the effect was assumed to follow from the center-of-mass motion of the pairs 
in the BEC side of the Feshbach resonance. However, in a 3D Fermi gas, also the 
RSA process produces the broadening effect. Furthermore, Fig. 4 in Ref. [8] shows 
asymmetric spectra for atoms of a given momentum. For a fixed momentum k the 
linear response theory predicts a symmetric spectrum centered at the single-particle 
resonance at — (Jt + -E* + The RSA diagram allows the transition with more 

energetic RF photons, producing a tail on the negative side of the single-particle peak. 

The present theory can be easily adapted for describing this experiment. Due to 
high temperatures, we assume that the initial state can be described by a normal finite 
temperature Fermi gas, implying that the quasiparticle energies in Eqs. (|12[) and (1131) 
can be replaced by single particle dispersion Ef. = |e& — fi\. The contributions from 
the two equations can be combined and we obtain in the limit of narrow linewidth 
(?7 — > 0) the transition probability 

\^k,p,q\W)\ = (^)2 + v >2 0{e k+q + € p - e k - e p - q - d), (20) 

where S(x) is the Dirac delta function and np(x) = 1/(1 + e^ x ~^>) is the Fermi-Dirac 
distribution at temperature T = l/fce/3, where fee is the Boltzmann constant. The 
lifetime broadening rf should in principle depend on the energy of the excitation, but 
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0.5 1 1.5 2 0.5 1 1.5 2 
Momentum fe/fe F Momentum k/k F 

Figure 5. The momentum resolved rf-spectra for the gapless system at zero 
temperature (left) and for T = 0.2 Tp (right). For low momenta k >C fcp, the 
spectrum follows the free particle dispersion. For high momenta k S> kp the 
spectrum follows a negative — 2 m dispersion. The transition region around 
fcp shifts towards higher momenta for increasing temperature in qualitative 
agreement with the experimental results. The white curves indicate the centers 
of gaussian fits to the individual spectra for different momenta k. These are to 
be compared with the experimental data in Refs. [8ll47|. 



here we will use the average lifetime obtained from the Galitskii equation [IB] 
rj 2(k F afId 3 p(l-^) 2 d(k F -\p\) 4 

E F tt J d 3 p6(k F - \p\) 5tt' (21) 

where we have replaced k F a = —2. This value is motivated by two-body scattering 
calculation which yields fc-dependent scattering amplitude [46]. It should be noticed 
that the qualitative results do not depend sensitively on the choice of this lifetime 
broadening. The momentum resolved spectra can be calculated by summing over 
momenta p and q in Eq. ([20]) and plotting the transition probabilities as functions of 
momentum k and energy e& — S. Doing this we obtain the spectra shown in Fig. [5j 
which can be directly compared with momentum resolved spectra in Refs. [8j 147) . 
Notice, however, that the latter experiment is performed at temperatures below the 
critical temperature, and it is somewhat surprising that in this case we obtain good 
match with the experimental spectra, having assumed a gapless state. This warrants 
a more detailed analysis, which will be a topic of further work. Despite the simplicity 
of the current model, the qualitative agreement is excellent. Clearly, normal state 
correlations are sufficient for producing the back-bending — dispersion, as also seen 
in Ref. [44] where a more complete T-matrix calculation was carried out for k ^> k F . 
The present model, however, can describe also the low momentum regime where the 
hole degrees of freedom become important. 

It is worth pointing out that possibly the most controversial assumption made in 
our theory, namely assuming that the atoms in hyperfme state |1) are described by 
bare propagators, was not needed in the above calculation of the momentum resolved 
spectrum for a normal state since the dressed and bare propagators are identical. 

In a recent experiment by the same group |48] the high momentum asymptote 



Resonant scattering effect in spectroscopies of interacting atomic gases 



15 



was studied also by measuring the expansion of the atom gas after abruptly switching 
off the interactions. This asymptote was found to match very well with the high 
momentum tail obtained from the rf-spectrum, in agreement with the theory of S. 
Tan [?TJ • It is worth noticing the similarity between transferring a particle into 
a noninteracting hyperhne state |3) and switching off of the interactions. While in 
the first the moment of the transition is not well defined (reflected as a time integral 
in the response), in the latter the moment when the interactions are switched off is 
determined by the external magnetic field. Despite this important difference, one can 
apply the RSA model also to the latter setting. The only contributing diagram is the 
Fig. [3] as the scattering must take place before the interactions are switched off. Now 
the corresponding integral formula is 

\ £ ^^ e -<^ +e >^^- e - g )( T -^^(. gi2C l, p c^ C2 , fc+9 c 1 , p _ 9 )|^(0))e^ t . (22) 

Evaluating this equation, and taking the limit T — > oo (one could also keep the time 
T finite but we want to compare the result with rf-spectra for which we considered 
the limit of infinitely long pulse), yields 

c\,p4,k c 2,k+ q c h p- q \ip(0)). (23) 



E l + £p + E l +q ~ ^ q + iV 



Now the probability of finding an atom in hyperfine state |2) with momentum k is 

9u J J 



^ El + el + El , „ - el „ + in x ' p 2 < 



Mm 



(24) 



,,,, *k + 
As above, in the discussion of the high energy tail of the rf-spectrum, we obtain 
constraints p — q^k + q^O and the energy change due to the scattering becomes 
E k~\~ e p~^ E k+ q ~~ e p- q ~ 2 £fc ' yielding the perceived occupation probability for momentum 
state k in the high momentum k — > oo limit 

512 2 1 (25) 



2e k + ir) k 4 ' 

This k~ 4 occupation probability produces the 8~ z / 2 asymptote in the rf-spectrum 
in linear response theory and this asymptotic behaviour is also exactly the same as 
obtained for rf-spectrum in RSA model. Thus the tail observed in Ref. [48] is also in 
agreement with the RSA model. However, we have made no assumptions regarding 
the occupation numbers for the initial state |^(0)). 



5.2. Bragg spectroscopy of a Fermi gas 

Another experiment used Bragg spectroscopy for studying a strongly interacting Fermi 
gas |12j . The experiment fitted very well with the theoretical mean- field predictions 
including both the heights and the widths of the spectral peaks. In the case of Bragg 
spectroscopy, the hyperfine state of the atom is not changed in the photon absorption 
but only the momentum state. Therefore, the hyperfine state |3) in the RSA diagrams 
corresponds to hyperfine state |2). Thus the final state interactions discussed in the 
context of rf-spectroscopy match the initial interactions 513 = 512 = g, and the 
diagrams interfere destructively. Thus one does not expect any contribution from 
the RSA effect to the Bragg spectroscopy in a Fermi gas. This is in good agreement 
with the experiment. This observation shows that the Bragg spectroscopy could be a 
good way for measuring single-particle properties in a Fermi gas since the RSA-type 
beyond mean-field effects are suppressed. 
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6. Comparison to experimental results: bosonic gases 

The RSA effect can be applied also to bosonic atoms. However, especially in the 
presence of a Bose-Einstein condensate, the effect of different diagrams is changed 
because atoms within the condensate do not feel the exchange interaction channel [55] . 
In the following, we will first consider a supcrfluid Bose gas and experiments performed 
on such systems, and later on discuss the Mott insulator state and the corresponding- 
experiments. 

6.1. RSA effect in BEC 

Let us consider a single-component Bose gas and let us assume, for the sake of 
simplicity, that all atoms are initially Bose-Einstein condensed into a zero-momentum 
state. Now the relevant RSA diagrams for Bragg and lattice modulation spectroscopies 
are the ones shown in Figs. [3] and [4] but since there is only one hyperfine component 
present, the hyperfine states |1), |2), and |3) are the same. The diagram in Fig. [2] 
does not contribute in the case of an ideal BEC because the diagram would require 
atoms initially in two different momentum states k and k + q and hence requires 
uncondensed atoms. Of course, such process would be allowed for non-ideal BEC and 
may be important in the case of strongly interacting BECs where quantum depletion 
is important or at high temperatures. 

As in the case of fermions, these bosonic RSA diagrams interfere coherently, but 
the only contribution comes from momenta k + q = p — q = 0. However, the diagram|4] 
experiences also the exchange interaction channel as the atom-atom scattering does 
not take place between two condensed atoms but between one condensed atom and 
one atom with momentum <7l- Hence the two RSA diagrams do not fully cancel each 
other. This partial cancellation of the RSA contribution produces interesting effects 
in the case of a superfluid BEC in an optical lattice. 

6.2. Lattice modulation spectroscopy of a Bose gas 

A superfluid Bose gas in a Id lattice was studied using lattice modulation 
spectroscopy [TS] and Bragg spectroscopy [H] [50] ■ In the experiments the atoms were 
transferred within the first lattice band [151 [50] as well as between the first lattice 
band and the second and third bands [14] . while changing the atom momentum due 
to the momentum carried by the photons. The transition to the second band was as 
predicted by the standard Bogoliubov (mean- field) theory, but the spectral peak of the 
transition to the third band and transitions within the first band were anomalously 
broad (see Fig. 3 in [TJ], Fig. 3a in [50J, and Fig.l in [H]). These effects were attributed 
to a loss of coherence [14] , quantum depletion of the condensate [15] , and to beyond 
mean-field correlation and thermal effects in a Id system [50 . 




The present theory suggests an alternative explanation. Since atoms in different 
lattice bands have reduced overlap, the interband interactions are reduced. In the 
harmonic potential approximation, atoms in the second band have an extra factor 
0.5 in the interactions with the first band (and in particular with the condensate), 
and the atoms in the third band have a factor 0.375. Including these factors, the 
bosonic RSA diagrams will cancel each other for the transition to the second band (as 
the contribution from the diagram @] is reduced by factor 0.5, cancelling the exchange 
interaction channel effect) but will yield only a factor 1 — 2 • 0.375 = 0.25 (contribution 
from diagram[4]is multiplied by 0.375) for the transition to the third band and a factor 





2J 




Figure 6. Bosonic RSA energy level scheme relevant for Bragg scattering and 
lattice modulation experiments |14| , |50| . and |15| . a) The atom-atom interaction 
U excites two particles (two dashed arrows) from the zero-momentum condensate. 
One of the atoms is then scattered by the photon S into the third band (the 
single dot-dashed line), b) In the Mott state, the atom-atom interaction U will 
yield a large energy shift due to the Mott gap, leading into a two-peak structure 
in the spectrum. The Bragg photon <5 provides an additional momentum kick 
(dot-dashed line) and needs to provide the energy for creation of the two Mott 
excitations. 



1 for the transitions within the first band. Thus the bosonic RSA effect does not affect 
the transitions to the second lattice band but it is present for the transitions to the 
first and the third bands. Similarly to the fermionic gases, the bosonic RSA effect 
will spread the spectral peaks. This is in qualitative agreement with the experimental 
observations. Furthermore, the width of the transition to the third band approaches 
the width of the third band and the width of the transitions within the first band is 
approximately twice the first band width, in agreement with present theory. Fig. [5] a) 
shows the schematic energy level diagram for the transition to the third band. 



6. 3. Bragg spectroscopy of a Bose gas 

The Bragg spectroscopy of a Bose gas has been studied also in a harmonic trap. 
The experiment in JILA in 2008 applied Bragg spectroscopy to study a strongly 
interacting BEC [Ml ESJ [51]. The width of the spectral peak showed anomalous 
broadening in the strongly interacting regime (see Fig. 3 in [E]). The bosonic RSA 
process explains the broadening by allowing transitions at higher energies as in the 
above discussion regarding the Id BEC experiments jT5] [T3J [50] and the momentum- 
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resolved rf-spectroscopy |S]. Moreover, the bosonic RSA process provides transitions 
also at lower energies, since non-condensed (at finite momentum k) atoms can be 
scattered into the condensate before absorbing the Bragg photon (in which case the 
atom-atom scattering before the photon absorption feels the exchange interaction 
channel but the opposite process does not). This reduces the kinetic energy of the 
atoms, allowing transitions at lower frequencies and thus further contributing to the 
width of the spectrum. These atom-atom scatterings down into the condensate are 
enhanced by the coherent state, and thus both RSA processes, scatterings up from the 
condensate and scatterings down into the condensate, are equally likely even when the 
quantum depletion is weak. Since the first process produces a high energy tail due to 
transitions to higher kinetic energy states, the latter process produces a similar low 
energy tail. Hence the total spectral lineshape remains symmetric, as also observed in 
the experiment [13) . 

Notice that a similar process as described by the RSA effect was crucial for the 
experiment, as the actual measured response was the total momentum of the atom 
gas [13] . The single particle excitations created by the Bragg scattering decayed due to 
scatterings with other atoms and the excess momentum yielded by the Bragg photon 
was thus transferred into center-of-mass momentum of the atom cloud. Thus the 
single-particle excitation was turned into a collective motion of the cloud. What our 
RSA effect now suggests, is that these atom-atom scatterings and the initial Bragg 
scattering are coherent in the sense that they are able to exchange energy as well as 
momentum. 

6.4- RSA effect in bosonic Mott insulator 

In the case of a Mott insulator state in a Bose gas in an optical lattice also the 
diagram [5] contributes as the condensate fraction vanishes and higher momentum 
states become populated in the ground state. Since there is no condensate present, 
all interactions experience the exchange interaction channel. However, the biggest 
difference to the interpretation of these diagrams as compared to supcrfluid gases 
is that the excitations created by the photon absorption do not couple to gapless 
excitations, as the gapless excitations in the Mott state are spin-flip excitations and 
these are not induced by atom-atom scatterings. Instead, the scattering of two atoms 
in the Mott insulator will necessarily create two single-particle excitations that are 
gapped by the Mott gap U, where U is the on-site interaction energy. The energy 
level scheme of the bosonic RSA process in a Mott insulator is shown in Fig. [6]b), 
showing how the process produces two Mott excitations. 

The Id Bose gases in the Mott insulator regime were studied using Bragg 
spectroscopy [50] and lattice modulation spectroscopy [15]. The measured spectra, 
such as in Figs.l f)-h) in [50] and Fig.l in [15], show two distinct peaks, the second 
resonance at twice the energy of the first resonance. The second peak was attributed to 
defects in the lattice, with atoms being transferred to doubly occupied sites [TS] . These 
observations fit in the picture drawn above as the additional peaks observed in |50] 
and [15] can be understood as a creation of two Mott excitations with energy cost 
equal to twice the Mott gap g. The process is made resonant by sufficiently energetic 
photon absorption, reflected as the second peak in the spectra. Notice that in the 
second order of atom-atom scatterings, one can also create three Mott excitations. 
Such processes should result in spectral peaks at three times the Mott gap. 
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6.5. Afterword 

To summarize, we analysed seven experiments in the light of our theory, and found 
out that there is a qualitative match. However, for the case of momentum resolved 
rf-spcctroscopy we performed quantitative analysis, obtaining good agreement with 
the experiment and the theory. As in the case of momentum resolved spectroscopy in 
Fermi gases, calculating quantitative predictions for also the other experiments using 
RSA-theory can be done. For example in the case of spectroscopies of superfluid 
BEC's in Refs. [H QH H ED], the starting point could be an ideal BEC and the 
interactions would then be considered only in the first-order perturbation theory by 
the present diagrammatic method. The approximation can be also easily improved by 
using the Bogoliubov theory, and the subsequent ground state, as the starting point. 
This is, indeed, an interesting topic for further work. 

7. Single-band Fermi Hubbard model 

Besides experiments, also exact numerical results provide a good platform for testing 
the theory. In this section we will study a Fermi gas in a one-dimensional optical 
lattice described by the single-band Fermi-Hubbard model. We will compare the 
predictions from the RSA theory to the exact rf-spectra obtained using time-evolved 
block decimation (TEBD) method. 

In the language used in these systems, the radio-frequency photon can be 
understood as creating a spin excitation and the RSA process describes then the 
decay of this excitation by scattering into lower energy spin excitation and creating a 
gapless charge excitation. The charge excitation thus plays the role of the Anderson- 
Bogoliubov phonon. We apply the time-evolved block decimation (TEBD) method as 
it allows us to access the exact radio-frequency spectrum of the system. The time- 
evolution of the system is determined by the Fermi-Hubbard Hamiltonian 

H(t) = -J 4a^+i<y+H.c. + uY,^2n l i+^e lSt/h Y,^ a c l2 +H.c (26) 

i,o-e{l,2,3} i i 

where U < is the attractive on-site interaction strength between atoms in states |1) 
and 1 2) and the density operator n ia — c\ a Ci ai where c^J destroys (creates) a particle 
of hyperfme state \a) in lattice site i. The state |3) is assumed to be noninteracting. 
Below we use units that have the hopping strength J = 1 . 

We use the TEBD algorithm for solving the ground state of the Fermi-Hubbard 
Hamiltonian in the absence of the radio-frequency field fi = and for empty hypcrfinc 
state 1 3). Starting from this ground state, we switch on the radio- frequency field and 
calculate the time evolution. Fig. [7] shows the fraction of atoms in state |3) after the 
radio-frequency pulse as a function of rf-detuning S. 

The BCS theory is known to fail in one-dimensional lattices. Even though the 
mean-field model does reproduce exact ground state energies faithfully, the behaviour 
of the pairing gap is phenomenally bad. In order to fix this, we do not solve 
the pairing gap A and the chemical potential /i from the BCS theory but rather 
derive them from the TEBD ground state following procedure in [52] . We use these 
values Atebd and /j.tebd in the BCS-type quasiparticle dispersion energies = 
\J [tk — /^tebd) 2 + A^ EBD , where the lattice dispersion = — 2Jcosfcforfc g [— tt, 7t], 
and in the Bogoliubov coefficients for the occupation numbers v\ — \ [\ — £fc ~^7 BBD ^ . 
The linear response rf-spectra, obtained using these TEBD derived parameters in 
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Figure 7. Exact radio-frequency spectra in one-dimensional lattice. The time- 
evolved block decimation method is used for calculating the exact rf-spectra for 
filling fractions (per spin state) a) 0.3 (1.1 J, —1.45 J), b) 0.5 (1.1 J, 0.0,7), 
and c) 0.7 (1.1 J, 1.45 J), with corresponding BCS fitting parameters shown as 
(A, fi), respectively. The on-site interaction strength is U = —5 J, rf-coupling 
strength Q = 0.025 J and the length of the pulse is T = 10 h/ J. Shown are also 
linear response (mean-field) spectra and spectra from the RSA model for both 
scatterings to higher and lower kinetic energy states (see text for discussion). 
The RSA transition probabilities are not properly normalized and hence the RSA 
spectra have been divided by 100 in order to better fit them in the picture. We 
have used 40 lattice sites in these simulations. 



the quasiparticle dispersion agree very well with parts of the observed TEBD 
spectra. However, the high energy tails of the spectra are missing from the linear 
response spectra, as the spectra are consistently too narrow. This shows that the 
BCS-type hole spectral function does not capture all of the relevant physics of the 
rf-photon absorption. The question is now whether this is a problem with the BCS- 
like description of the ground state (that is, how well does the BCS-like quasiparticle 
dispersion work even when used with the TEBD derived parameters) or with the 
dynamics described by the linear response theory. With this caveat in mind, let us 
consider how would the RSA process affect the picture. 

Fig.[7]shows also the RSA spectra from the Eqs.[T2land[T51 This contribution has 
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been separated into two parts, one corresponding to atom-atom scatterings to higher 
kinetic energy states AK kpq = E\ +q + e p — E\ — e p _ q > and one for scatterings to 
lower kinetic energy states AKf. pq < 0. While the former process can be understood 
as a phonon excitation, the latter process may be unphysical and possibly an artifact 
from our approximative approach. We will neglect these scatterings to lower energy 
states for a moment and discuss them later on. 

The combination of the standard linear response spectrum and the RSA spectrum 
seems to describe rather well all qualitative features of the exact TEBD spectrum: 
Eq. 12.11 reproduces the main peak of the exact spectrum and Eq. [T^] provides the 
correct width for the spectrum (neglecting the AK < contribution). Indeed, the 
8 J width of the spectrum can be observed in all TEBD spectra at different filling 
fractions and different interaction strengths (although resolution becomes difficult for 
weak interactions — U < 4 J). This width in the one-dimensional lattice is also a 
general property of Eqs. (fT2"j) and (HUl) . The low energy end of the spectrum is provided 
by low momentum \2,k — 0)-atoms that scatter from \l,p = fcp)-atoms at the Fermi 
surface but the change in momentum is small (q — 0). Effectively it corresponds to 
no scattering at all, and the resonant energy for the transition is 5\ ow — £q — \i + Eq . 
The high energy end of the spectrum is provided by |2, k = it) -atoms at the edge of 
the lattice band scattering from \l,p = 0) -atoms in the bottom of the Fermi sea with 
momentum exchange q — — it, yielding resonant energy <5hi g h = + Eq + e\ — e\ — [i. 
Thus the total width of the spectrum equals Jhigh — <>iow = e„ + — to — e o = 8 J, 
independent of the chemical potential /i (and hence filling fraction) and the excitation 
gap A. This is in excellent, quantitative agreement with the exact TEBD spectra. 
Notice also that this implies that the width of the spectrum is actually insensitive 
to the way how the A and [i are determined, and even the values from the BCS 
theory would yield the same total width; only the positions of the peaks would be 
affected. And finally, the quasiparticle dispersions vanish from the calculation of 
the total width, and hence even if one would relax the assumption of the BCS-type 
quasiparticle dispersion E^ somewhat, the result would remain the same. 

At high filling fractions our model predicts a low energy tail extending beyond the 
8 J width in disagreement with the TEBD spectra. The tail is caused by scatterings 
to lower energy states AK kpq < 0. Since the |l)-atoms are forbidden from scattering 
to lower kinetic energy states due to the presence of a well-formed Fermi sphere, the 
lowering of the energy is possible only if the increase in the kinetic energy of the 
|l)-atom is lower than the decrease in the quasiparticle energy of the |2)-atom. Thus 
we have condition 

> AK kpq = (E k+q - E k ) + (e p _ g - e p ). (27) 
In the long scattering wavelength limit q — > this yields 
dE k de p 

-^<-^U=fe- ( 28 ) 

where the right-hand side must be evaluated close to the Fermi surface because of the 
well-formed Fermi sphere for |l)-atoms. The left hand side yields e ^T M ^ , where the 
prefactor (e k — \i)jE k has magnitude less than 1. However, since the quasiparticles 
can populate any momentum state k, it is possible that Eq. (|28|) is satisfied for some 
k as long as the single particle dispersion is not too rapidly increasing function at 
the Fermi momentum fcp- In particular, at high fillings the single-particle dispersion 
e k becomes concave at the Fermi surface and the scatterings to lower kinetic energy 
states AK < becomes possible. At half filling, the single-particle dispersion has 
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the highest slope and no AK < scatterings are possible. At low fillings there are 
some scattering channels as well but the contribution is very small as compared to the 
scatterings to higher energy states AK > (typically at least one order of magnitude 
lower) . 

8. Discussion 

In conclusion, we have formulated a first order scattering theory which adds to the 
understanding of interdisciplinary spectroscopies. Despite being of the first order 
in the photon coupling, the theory goes beyond standard mean-field linear response 
theories by incorporating an approximative coupling to phonon-like modes. Similar 
phenomena have been extensively studied in nuclear physics since 1960's [35] . and it 
would be interesting to compare the present approach with the models used in nuclear 
physics. The resonant scattering process in ultracold atom gases becomes increasingly 
important with stronger interactions and its effect can be seen in a wide range of 
experiments - here we have analyzed seven of them. 
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